Physics of Fluids 
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To study subregions of a turbulence velocity field, a long record of velocity data of grid turbulence 
is divided into smaller segments. For each segment, we calculate statistics such as the mean rate 
of energy dissipation and the mean energy at each scale. Their values significantly fiuctuate, in 
lognormal distributions at least as a good approximation. Each segment is not under equilibrium 
between the mean rate of energy dissipation and the mean rate of energy transfer that determines 
the mean energy. These two rates still correlate among segments when their length exceeds the 
correlation length. Also between the mean rate of energy dissipation and the mean total energy, 
there is a correlation characterized by the Reynolds number for the whole record, implying that the 
large-scale flow affects each of the segments. 



I. INTRODUCTION 

For locally isotropic turbulence, Kolmogorov^ consid- 
ered that small-scale statistics are uniquely determined 
by the kinematic viscosity v and the mean rate of energy 
dissipation (e). The Kolmogorovvelocity uk = (I'is))^^^ 
and the Kolmogorov length -q = (z/'^/(e))^/^ determine 
the statistics of velocity increment (5u,- = u{x -\-r) — u{x) 
at scale r as 

^ = F„f-) for n = 2,3,4,.... (1) 

Here (•) denotes an average over position x, and Fn is 
a universal function. The universality is known to hold 
well. While (^u") at each r is different in different veloc- 
ity fields, (e) and hence and -q are accordingly differ- 
ent. That is, (e) is in equilibrium with the mean rate of 
energy transfer that determines {Su^). 

However, the universality of small-scale statistics 
might not be exact. To argue against the exact univer- 
sality. Landau^ pointed out that the local rate of energy 
dissipation e fluctuates over large scales. This fluctu- 
ation is not universal and is always significant i^i^ In 
fact, the large-scale flow or the configuration for tur- 
bulence production appears to affect some small-scale 
statistics,^i^ii^i^ii2, 

Obukho-^^ discussed that Kolmogorov's theory^ still 
holds in an ensemble of "pure" subregions where e is 
constant at a certain value. Then, the e value represents 
the rate of energy transfer averaged over those subre- 
gions. For the whole region, small-scale statistics reflect 
the large-scale flow through the large-scale fluctuation of 
the £ value. The idea that turbulence consists of some 
elementary subregions is of interest even nowi^ We study 



''Electronic address: |hmourl@mri-jma.go.jp| 

''Also at Meteorological and Environmental Sensing Technology, 

Inc., Nanpeidai, Ami 300-0312, Japan 

'^Electronic address: mtakaoka@maiI. doshisha. ac.jp| 



statistics among subregions in terms of the effect of large 
scales on small scales, by using a long record of velocity 
data obtained in grid turbulence. 



II. EXPERIMENT 

The experiment was done in a wind tunnel of the Mete- 
orological Research Institute. Its test section had the size 
of 18, 3, and 2 m in the strcamwise, spanwise, and floor- 
normal directions. We placed a grid across the entrance 
to the test section. The grid consisted of two layers of 
uniformly spaced rods, with axes in the two layers at right 
angles. The cross section of the rods was 0.04 x 0.04 m^. 
The separation of the axes of adjacent rods was 0.20 m. 

On the tunnel axis at 4 m downstream of the grid, 
we simultaneously measured the streamwise {U -\- u) and 
spanwise {v) velocities. Here U is the average while u{t) 
and v{t) are fluctuations as a function of time t. We 
used a hot-wire anemometer with a crossed-wire probe. 
The wires were made of platinum-plated tungsten, 5 /im 
in diameter, 1.25 mm in sensing length, 1mm in sep- 
aration, oriented at ±45° to the streamwise direction, 
and 280 °C in temperature. The signal was linearized, 
low-pass filtered at 35 kHz, and then digitally sampled at 
fs — 70 kHz. We obtained as long as 4 x 10® data. 

The calibration coefficient, with which the flow veloc- 
ity is proportional to the anemometer signal, depends 
on the condition of the hot wires and thereby varied 
slowly in time. We determine the coefficient so as to 
have U = 21.16ms~^ for each segment with 4 x 10® data. 
Within each segment, the coefficient varied by ±0.4% 
at most. Also varied slowly in time the flow temper- 
ature and hence the kinematic viscosity v. We adopt 
V — 1.42 X 10^^ m^ s~^ based on the mean flow temper- 
ature, 11.8°C. The temperature variation, ±1.2 °C, cor- 
responds to the V variation of ±0.7%. These variations 
are small and ignored here. 

Taylor's frozen-eddy hypothesis, i.e., x — —Ut, is used 
to obtain u{x) and v{x) from u{t) and v{t). This hypoth- 
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TABLE I: Turbulence parameters: mean energy dissipa- 
tion rate (e), rms velocity fluctuations {u^}^^^ and {v^)^^"^ , 
Kolmogorov velocity uk, flatness factors (u'*)/(u^)^ and 
{v^)/{v^y, correlation lengths Lu, Ly, and L^, Taylor mi- 
croscale A, Kolmogorov length rj, and microscale Reynolds 
number Re^. 



Quantity 



Value 



(He)) 



1/4 



= /o°°(«(^ + '-)«(a;)><ir/{u2) 
Lt, — Jq°° {v{x + r)v{x))dr / (f ^) 

= f„°°{e(x + r)e{x) - {ef)dr/{s 
X = [2{v^)/{(a^vf)]^/^ 



7.98 m^s"-" 
l.lOms-i 

i.oems-"^ 

O.lOSms^' 
3.02 
3.00 
17.9 cm 
4.69 cm 
0.469 cm 
0.548 cm 
0.0138 cm 
409 



esis requires a small value of {u^y^^/U. The value in our 
experiment, 0.05, is small enough. Since u{t) and v(t) 
are stationary, u{x) and v{x) are homogeneous, although 
grid turbulence decays along the streamwise direction in 
the wind tunnel. We are mostly interested in scales up 
to about the typical scale for energy-containing eddies, 
which is much less than the tunnel size. Over such scales, 
fluctuations of u{x) and v{x) correspond to spatial fluc- 
tuations that were actually present in the wind tunneli^ 
Those over the larger scales do not. They have to be in- 
terpreted as fluctuations over long timescales described 
in terms of large length scales 

Turbulence parameters are listed in Table IH Here and 
hereafter, (■) is used to denote an average over the whole 
record. The derivative was obtained as dxV = [8v{x + 
Sx) - 8v{x - 6x) - v{x + 2Sx) + v{x - 2Sx)]/12Sx with 
Sx — U/fs- The local rate of energy dissipation was 
obtained as e = 15i'{dxv)^ /2 instead of usual 15i'{dxu)^ , 
in order to avoid possible spurious correlations with 6ur 
over small r for analyses in the next section. 

Figure [1] shows {SuD/u'^, the u, v, and e correlations, 
and also the correlation lengths L^, L^, and L^. We see 
the inertial range, albeit narrow, where {Su'^) roughly 
scales with r^/"^ . The u and v correlations are significant 
up to r ~ lO'^ry, which corresponds to the scale of largest 
eddies. The correlation length L„ corresponds to the typ- 
ical scale for energy-containing eddies. Since e belongs 
to small scales, its correlation decays quickly. 



Obukhov.-'^^ For each segment, we have statistics such as 

idxv)l{x) 
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v(x' -f r)v{x')dx' 
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Here x is the center of the segment, and r < R. The 
mean rate of energy dissipation is £_r, — 15v{dxv)'j^/2, 
which yields the Kolmogorov velocity uk.r — {lySRY^^ 
and the Kolmogorov length rj^ — (y^ j e b)^ ^ . We also 
have the mean total energy, v\ — v(x -\- r)v(x)R for 
r — 0, and the microscale Reynolds number, RcA.i? = 
2^/^w|j/j/[(92.z;)^]-^/^. The mean rate of energy trans- 
fer, however, is not available from our experimental data. 
Fig. [Ija) shows an example of Suf, j^/u^ ^, which differs 
from {5ul)/ul^. 



A. Distribution of fluctuation 

Over a range of i?, we study statistics of Su^ ^/ mI- ^ 
among segments. The scale r is fixed at IOt^^^ in the 
dissipation range and IQQrjR in the inertial range. Since 
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III. RESULTS AND DISCUSSION 

The data record is now divided into segments with 
length R. They correspond to subregions considered by 



FIG. 1: (a) {Suf.)/u^ as a function of r/r] and Su^ j^/u'^ ji for 
a segment with R = W^ri as a function of r/rjR. We indicate 
Z/u, Lv^ L^, and the inertial range. 

{v{x + r)v{x))/{v'^), and {e{x + r)e{x) - {ef)/{e^ - {ef) as 
a function of r/ry. 
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6ur is available only at discrete scales r that are multiples 
of the sampling interval [///«, Suj^^^ j^ and Suj^Q^^ j^ 
are obtained through interpolation by incorporating the 
fluctuation of 777? among segments. 

Figure ^a) shows the standard deviation. The fluc- 
tuation of Su^ r/uk _r at a fixed r/ijn is significant even 
when R is large. In individual segments, Kolmogorov's 
theoryi does not hold. The mean rate of energy transfer 
that determines Su^ ^ is not in equilibrium with the mean 
rate of energy dissipation e r that determines ^ and 
77fl. The degree of this nonequilibrium fluctuates among 
segments and thereby induces the observed fluctuation of 

The mean rate of energy transfer also fluctuates among 
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FIG. 2: Statistics of 5uior,H,fl/'iK,fl (filled triangles down), 
'^■"ioo?7h,«/'"k,h (filled triangles up), (open circles), vj^ 
(open squares), and ReA,ij (open diamonds) as a function 
of R/rj. (a) Standard deviation normalized by the average. 
The dotted line indicates the R~^^^ scaling, (b) Skewness 
factor, (c) Flatness factor. Also shown are the skewness and 
flatness factors of ^n{5uiQ^^ [i/u^ji) (open triangle down), 
ln((5uiooT;n,fl/'^K,R) (open triangles up), Inefl (filled circles), 
\nv^ (filled squares), and InRcA.ij (filled diamonds). 
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FIG. 3: (a) Probability densities of Xr = 5MioO)7R,fl/*^K,H 
(triangles up) and ln(5ufoO)7H,ii/^K.ii) (triangles down) at 
R = 10^77 (open triangles) and 10^»7 (filled triangles) as 
a function of {Xr - {Xr))/{{Xr - {Xr))^^^- Those of 
Xr — SuIq^^^ r/uj^^r are shifted by a factor of 10^. The 
dotted lines denote the Gaussian distribution, (b) Same as 
(a) but for Xr = er (triangles up) and Insij (triangles down). 



scales in each segment, which does not necessarily have 



the inertial-range scaling Sw^ ^ 



p2/3 



[Fig. ma) 



If each segment had this scaling, the fluctuation of 
Sulfj/u^j^ among segments at a fixed r/rjn in the in- 
ertial range would correspond to the fluctuation of the 
Kolmogorov constant Su^ ^/ {rea)^^^. 

Figures m^b) and[2l|c) show the skewness and flatness 
factors of ln((5u^ r/^k r) (open triangles). They are close 
to the Gaussian values of and 3. Thus, at least as a 
good approximation, the distribution of Su'^j^/u^j^ is 
lognormal. This is also the case in er (filled circles) and 
at i? > lO'^ry ~ L„ in vj^ and ReA,fl (filled squares and 
diamonds), while the mean rate of energy transfer should 
not have a lognormal distribution because it changes its 
sign. Examples of the probability density functions are 
shown in Fig. [3l The lognormal distribution of er was 
discussed as a tentative model by Obukhov.^^ 

A lognormal distribution stems from some multiplica- 
tive stochastic process, e.g., a product of many inde- 
pendent stochastic variables with similar variances. To 
its logarithm, if not too far from the average, the cen- 
tral limit theorem applies. For the lognormal distribu- 
tions observed here, the process is related with the en- 
ergy transfer. While the mean energy transfer is to a 
smaller scale and is significant between scales in the in- 
ertial range alone, the local energy transfer is either to 
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FIG. 4: (a) (<5u?,fl/uic,fl)e at = 10^ 7) as a function of 
r/rjn (solid lines) and {5u^)/u^ as a function of r/77 (dotted 
line). The inset shows the dependence of {5u\q^^ji/vj^j))^ 
on {eR)^/{e). The horizontal dotted line indicates the value 
of {Sulori)/uK- (b) Same as (a) but at J? = 10^ r;. 



FIG. 5: (a) (v\)i,/{v'^) as a function of (e_R)j/(e) at 
R — 10^77 (triangles), 10^77 (circles), and 10*7; (squares). 
Along the dotted lines, the microscale Reynolds number 
1^^^^ {v\)s / {^{eR)^)^^^ is constant, (b) {R.e\,R)e as a func- 
tion of {eR)e/{£)- The symbols are the same as in (a). 



a smaller or larger scale and is significant between all 
scales i^iiSiii^ Any scale is thereby affected by itself and 
by many other scales. They involve large scales because 
the lognormal distributions are observed up to large R. 
There is no dominant effect from a few specific scales, in 
order for the central limit theorem to be applicable. 

At R > 10^77 ~ 10^i„, there are alternative 
featuresi^ii^ The standard deviations scale with R^^^^ 
[Fig. Uta)]. The skewness and flatness factors of 
'^'"r _r/"k fl' ^-R' '"r^ ^^'^ R.eA,_R are close to the Gaus- 
sian values [Figs. [2jb) and[2{c): filled triangles, open 
circles, open squares, and open diamonds; see also Fig. 

Their distributions are regarded as Gaussian rather 
than lognormal, although this has to be confirmed in fu- 
ture using high-order moments or probability densities at 
the tails for the larger number of segments. The R^^^^ 
scaling and Gaussian distribution are typical of fluctua- 
tions in thermodynamics and statistical mechanics, for 
which no correlation is significant at scales of interest as 
in our case at r > 10^77 [Fig. mb)]. 



B. Correlation between fluctuations 

The fiuctuations among segments at i? > in have in- 
teresting correlations. Since these correlations are weak, 
they are extracted by following Obukhov,^^=. i.e., by av- 



eraging over segments with similar er. Specifically, we 
use conditional averages, denoted by (•)£, for ranges of 
er separated at (e)/4, {e)/2, (e), 2(e), and 4(e). 

Figure [3] shows (<^Ur fl/'"K i?)^ ^ ^ function of r/rjR. 
When R =. 103,7 ~ L„' [Fig.'Hb)], (SulJulj^), is in- 
dependent of {sfj)^. The former varies only by a factor 
of 1.5 while the latter varies by a factor of 20. When 
R = 102,7 ~ 10-iL„ [Fig. Ha)], {SulJulj^}, is not 
independent of {eji)^. 

The implication of the above result is that the mean 
rate of energy transfer that determines Su^ ^ correlates 
with the mean rate of energy dissipation er that deter- 
mines and rjR among segments with R^ Lu- Here 
Lu is the typical scale for energy-containing eddies. Most 
of the energy of such an eddy is transferred through scales 
and dissipated within its own volume. Thus, each energy- 
containing eddy tends toward equilibrium between the 
mean rates of energy transfer and dissipation. This ten- 
dency does not exist at i? < L^-^^ Within an energy- 
containing eddy, the spatial distribution of c is not ho- 
mogeneous. In fact, the e correlation is significant at 
r<Lu [Fig. IHb)]. 

Therefore, in order for statistics such as Su^j^/u^j^ 
and ((Jm^ _r/"k have physical meanings, the mini- 

mum segment length is about L„. We are to discuss that 
such segments individually reflect the large-scale flow. 

Figure [Ha) shows {vr)s/{v'^) as a function of 
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{eii)e/{e). The former varies with the latter. When 
R > 10^77 ~ Lu, the microscale Reynolds number 
^ is almost constant at the value for the 

whole record, Re\ = 409. Thus, segments obey a correla- 
tion between vj^ and en characterized by the ReA value. 
Since Rca is determined by the large-scale flow, it follows 
that the large-scale flow affects each of the segments. A 
consistent result is obtained for (ReA.fi)^ in Fig. [^b). 

The above tendency toward a constant microscale 
Reynolds number originates in energy-containing scales. 
In general, through an empirical relation (e) oc 
{v'^)^/'^ /Ly, Rca is related to the Reynolds number for 
energy-containing scales as Rca oc {{v'^y^'^Ly/vY^^. 
Then, Rca is constant if (x (v^)^^^^. Fig. [6] shows 
{v{x + r)v{x)i^)s/ {vj^)i;, which extends to larger scales 
for smaller (vj^)^. The correlation length {v{x + 
'i^)v{x)ii)edr/{v'j^)i. should be accordingly larger. 

The process for each segment to reflect the large-scale 
flow could be related to the energy transfer. As noted 
before, the energy transfer couples all scales. The energy 
transfer itself is affected by the large-scale flow. This 
is because energy is transferred between two scales via 
an interaction with some other scale. When the inter- 
action occurs with a large scale, the energy transfer is 
strong . 



C. Unimportance of e fluctuation 

Obukho\J^ discussed that, through the fluctuation of 
the en value, the large-scale flow affects small-scale statis- 
tics for the whole region such as {Su")/u^. The reason 
is that is obtained at a fixed r, regardless of the 

fluctuations of uj^ ^ and rjn induced by the fluctuation of 
en. We are to discuss that the fluctuation of the en value 
is not important so far as i? L^. This condition on R 
is required for the en value to correlate with and thus 
statistically represent the mean rate of energy transfer 
that determines Su^ n- 

Figure [7] compares ((Su")/-Uk to (Sul^j^/u^ j^), i.e., av- 
erage of Su^ j^/u^ ^ at each r/rjn over all segments, with 
R = lO^r; ~ Lu. Even for n = 6 [Fig. [TJb)], they are 
not distinguishable [see also Fig. IUb)]. Although en 
and hence UK,n f^nd rjn fluctuate among segments, these 
fluctuations are not large enough for {Su^)/u^ to differ 
from i^u" j^J j^) . This conclusion is general. In vari- 
ous flows,^!^ including an atmospheric boundary layer 
at Rca — 9000, the standard deviation of en/i^R) at 
R ~ Lu is close to the value obtained here [Fig. HUa)]. 
Hence, through the fluctuation of the en value, the large- 
scale flow does not affect small-scale statistics for the 
whole record. It was suggested that the large-scale flow 
does affect {6u^)/u^, even in the scahng exponentsi^'^i'i^ 
If this is the case, the effect is already inherent in the 
individual segments. They are unlikely to be "pureii" or 
elementary. 
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(solid line) and {Su")/u^ as a function of r/rj (dotted line), 
(a) n = 2. (b) n = 6. 



IV. CONCLUDING REMARKS 

Using segments of a long record of velocity data ob- 
tained in grid turbulence, we have studied fluctuations 
of statistics such as ^u^j^/u^j^, en, and The fluc- 
tuations are signiflcant and have lognormal distributions 
at least as a good approximation (Figs. [5] and [3]). In 
each segment, the mean rate of energy transfer that 
determines Su^ ^ is not in equilibrium with the mean 
rate of energy dissipation en that determines ^ and 
rjn- These two rates still correlate among segments with 
R ^ Lu (Fig. 3]) , which tend toward equilibrium between 
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the two rates. Also between er and t;^, there is a corre- 
lation characterized by Re^ for the whole record (Fig. O. 
Thus, the large-scale flow affects each of the segments. 

The observed fluctuations depend on L„ and Rca, 
which in turn depend on the configuration for turbu- 
lence production, e.g., boundaries such as the grid used 
in our experiment. Nevertheless, the significance of those 
fluctuations implies that they have been developed in 
turbulence itself. Their lognormal distributions are ex- 
plained by a multiplicative stochastic process in turbu- 
lence, which is related with the energy transfer among 
scales. The correlations among the fluctuations are also 
explained in terms of the energy transfer. 

Previous studies^'^'^iiSiiSi'i^ suggested that the large- 
scale flow affects some small-scale statistics, although this 
has to be conflrmed at higher Rba where large and small 
scales are more distinct. If the effect really exists, it is 
inherent individually in the segments. 

Our study was motivated by Obukhov's discussion. 
It implies the presence of equilibrium between the mean 
rates of energy transfer and dissipation in an ensemble of 
segments with similar values of er. This is the case at 
R }^ Lu (Fig. 2]).^° Also as discussed by Obukhov, the 
distribution of er is lognormal at least as a good approxi- 
mation (Figs. [2]and[3]). However, although Obukhov dis- 
cussed that the large-scale flow affects small-scale statis- 
tics through the fluctuation of the er value, this is not 
the case (Fig. [T]). 

The lognormal distributions observed here have to 
be distinguished from those proposed by Kolmogorovi^i 
While he was interested in small-scale intermittency and 
studied Er and (5u" at small r to obtain their scaling 
laws, we are interested in large-scale fluctuations and 
have studied er and Ju" ^ at small r but at large R. The 



scaling laws of er and are not necessary. In addi- 

tion, the lognormality has been attributed to a different 
process. Hence, our study is not necessarily concerned 
with the well-known problems of Kolmogorov's lognor- 
mal model, e.g., violation of Novikov's inequality^^ for 
scaling exponents. Still exists a possibility that small- 
scale intermittency is affected by large-scale fluctuations. 
The study of this possibility is desirable. 

There were no studies of statistics among segments 
with large R. Hence, we have focused on grid turbu- 
lence, which is simple and thus serves as a standard. 
For flows other than grid turbulence, the fluctuations of 
statistics among segments are expected to be signiflcant 
as well. In fact, regardless of the flow conflguration and 
the Reynolds number, the large-scale fluctuation of er 
is signiflcant.'^''''^ Those fluctuations are also expected to 
have lognormal distributions and mutual correlations as 
observed here because they are due to the energy trans- 
fer in turbulence itself. However, grid turbulence is free 
from shear. It was previously found that Sur correlates 
with u in shear flows such as a boundary layer but not in 
shear-free flows. ^'^'^ The fluctuations of statistics among 
segments might be somewhat different in a shear flow. 
To this and other flow conflgurations, it is desirable to 
apply our approach. 
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